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Abstract. We present our recent results on the Prandtl-Meyer reflection for 
supersonic potential flow past a solid ramp. When a steady supersonic flow passes 
a solid ramp, there are two possible configurations: the weak shock solution 
and the strong shock solution. Elling-Liu's theorem (2008) indicates that the 
steady supersonic weak shock solution can be regarded as a long-time asymptotics 
of an unsteady flow for a class of physical parameters determined by certain 
assumptions for potential flow. In this paper we discuss our recent progress in 
removing these assumptions and establishing the stability theorem for steady 
supersonic weak shock solutions as the long-time asymptotics of unsteady flows 
for all the physical parameters for potential flow. We apply new mathematical 
techniques developed in our recent work to obtain monotonicity properties and 
uniform apriori estimates for weak solutions, which allow us to employ the Leray- 
Schauder degree argument to complete the theory for the general case. 



1. Introduction 

We are concerned with unsteady global solutions for supersonic flow past a solid 
ramp, which may be also regarded as portraying the symmetric gas flow impinging 
onto a solid wedge (by symmetry). When a steady supersonic flow past a solid 
ramp whose slope is less than a critical slope, Prandtl employed the shock polar 
analysis to show that there are two possible configurations: the weak shock reflec- 
tion with supersonic downstream flow and the strong shock reflection with subsonic 
downstream flow, which both satisfy the entropy conditions, provided that we do 
not give additional conditions at downstream; see Busemann [3], Courant-Friedrichs 
[ID], Meyer [15], and Prandtl [IB] . 

The fundamental question of whether one or both of the strong and the weak 
shocks are physically admissible has been vigorously debated over the past seventy 
years, but has not yet been settled in a definite manner (cf. Courant-Friedrichs 
|10| . Dafermos [11], and Serre [II])- On the basis of experimental and numerical 
evidence, there are strong indications that it is the Prandtl-Meyer weak reflection 
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solution that is physically admissible. One plausible approach is to single out the 
strong shock reflection by the consideration of stability, the stable ones are physical. 
It has been shown in the steady regime that the Prandtl-Meyer weak reflection 
is not only structurally stable (cf. Chen-Zhang-Zhu [8]), but also L 1 -stable with 
respect to steady small perturbation of both the ramp slope and the incoming steady 
upstream flow (cf. Chen-Li [7]), while the strong reflection is also structurally stable 
for a large spectrum of physical parameters (cf. Chen- Fang [9]). The first rigorous 
unsteady analysis of the steady supersonic weak shock solution as the long-time 
behavior of an unsteady flow is due to Elling-Liu in [13] in which they succeeded 
to establish a stability theorem for a class of physical parameters determined by 
certain assumptions for potential flow (see §3.1). 

The purpose of this work is to remove the assumptions in Elling-Liu's theorem [13] 
and establish the stability theorem for the steady supersonic weak shock solution as 
the long-time asymptotics of an unsteady flow for all the admissible physical param- 
eters (without additional conditions) for potential flow. To complete this theorem, 
we apply new mathematical techniques developed in Chen-Feldman [6] to obtain 
uniform apriori estimates for weak solutions. We first establish the monotonicity 
property and its consequence of weak solutions (see §4.1). Then we make various 
uniform apriori estimates of weak solutions for two cases: the C 2 '°-estimates away 
from the sonic circles where the governing equation is uniformly elliptic (see §4.3.1) 
and the weighted C 2,Q -estimates near the sonic circles where the ellipticity degen- 
erates (see §4.3.2). These careful estimates allow us to employ the Leray-Schauder 
degree argument to establish the complete theory (see §4.4). 

2. Mathematical Formulation for the Problem 

The unsteady potential flow is governed by the conservation law of mass and 
Bernoulli's law: 

dtp + V x • (pV x $) = 0, (2.1) 

dt$ + \\V x <S>\ 2 + i(p) = B (2.2) 

for the density p and the velocity potential where the Bernoulli constant B is 
determined by the incoming flow and/or boundary conditions, and i(p) satisfies the 
relation 

i'( P ) = m = fipi 
p p 

with c(p) being the sound speed, and p is the pressure that is a function of the 
density p. For an ideal polytropic gas, the pressure p and the sound speed c are 
given by 

p(p) = Kp' 1 , c 2 {p) = Kjp' 1 ' 1 
for constants 7 > 1 and k > 0. Without loss of generality, we choose k = I/7 to 
have 

i(p) = ^—=r^, c\p)=p^\ (2.3) 
7-1 

This can be achieved by the following scaling: 

(x.,t,B) — > (ax, a 2 t, a~ 2 B), a 2 = kj. 



Taking the limit 7 — > 1+, we can also consider the case of the isothermal flow 
(7 = 1). For isothermal flow, f|2.3|) implies 

i(p) = In p, c 2 {p) = 1. 

Our goal is to find a solution (p, $) to system (|2.ip - (|2.2p when, at t = 0, a uniform 
flow in R\ := {x x € R,x 2 > 0} with (p, V x $) = (p ooj^oojO) is heading to a solid 
ramp: 

W := {x = (x\, X2) : < X2 < x\ tan 9 W , x\ > 0}. 



Problem 1 (Initial-Boundary Value Problem). Seek a solution of system (|2.ip - (|2,2p 

with B = + and i/ie initial condition at t = 0: 

(p,$)|t=o = (Poo,^ooXi) /or (xi,x 2 ) €R+\W, (2.4) 
and mi/i the slip boundary condition along the wedge boundary dW: 

Vx$ • v\dWn{x 2 >0} = 0, (2.5) 
where v is the exterior unit normal to dW . 

Notice that the initial-boundary value problem (f2TT]) - (p^]) with ([2T4]) - ([23]) is 
invariant under the scaling: 

(x, t) — > (ax, at), (p, <£) — > (p, — ) for 

a 

Thus, we seek self-similar solutions in the form of 

p(x, t) = p{£, rj), $(x, t) = t<f>{£, rj) for (£, 7?) = ~. 

Then the pseudo-potential function tp = <f> — \ {^ 2 + rj 2 ) satisfies the following Euler 
equations for self-similar solutions: 

div(pDp) + 2p = 0, (2.6) 

£ — + -\Dip\ 2 + p) = B, (2.7) 

7 — 1 I 

where the divergence div and gradient D are with respect to (£, n). From this, we 
obtain an equation for the pseudo-potential function (/?(£, rj) as follows: 

div(p(|^| 2 , <p)Dtp) + 2p(\Dip\ 2 , ip) = (2.8) 

for 

p{\Dp\ 2 ^) = (B - (7 - l)(\\D V \ 2 + v))^\ (2.9) 
where we set Bq := (7 — 1)S + 1. Then we have 

c 2 (|ZV| 2 , ¥>) = 5 - (7 - ^(^l^l 2 + ( 2 - 10 ) 

Equation (|2.8p is an equation of mixed elliptic-hyperbolic type. It is elliptic if and 
only if 

\Dp\ <c(\D^p\ 2 ,p). (2.11) 
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If p is a constant, then, by (|2.8j) and (|2.9p . the corresponding pseudo-potential ip 
is in the form of 

»7) = " g + ^) + < + ^ + k 

for constants u, u, and k. As the incoming flow has the constant velocity (uoo,0), 
the corresponding pseudo-potential tpoc has the expression of 

Poo = —it 2 + V 2 ) + Uoct + koo (2.12) 

for a constant koo. 

Without loss of generality by scaling, we fix p^ = 1 and = (7 "_°^ )/2 = 

Pec 

(Mach number of the state at infinity), and then (|2.7p becomes 

Then Problem [1] can be reformulated as the following boundary value problem in 
the self-similar coordinates (C,v)- The domain in the self-similar coordinates (£,7/) 
corresponding to {(x, t) : x £ \ W, i > 0} is 

A:=Rl\{^, V ):r]<Ct a n9 w ,C>0}. 



Problem 2 (Boundary Value Problem) . Seek a solution ip of equation (12. 8p in the 
self-similar domain A to^/i the slip boundary condition: 

D(p-v w = on T wedge = {(£,77) : n = £ tan 9 W ,£ > 0}, (2.14) 

where v w is the exterior unit normal to the boundary of the wedge ^ wedge- 

In particular, we seek a weak solution of Problem[2]such that it contains a straight 
weak oblique shock attached to the tip of the wedge, and the oblique shock is 
connected to a normal shock through a curved shock as shown in Figured! 




Figure 1. Weak shock solutions in the self-similar coordinates 

A shock is a curve across which Dp is discontinuous. If f2 + and J7 _ (:= \ S7+) 
are two nonempty open subsets of O C M 2 , and S := d£l + n is a C 1 -curve where 
Dp has a jump, then p £ W/o' c ^ C ,1 (0 ± U 5) fl C 2 (J7 =t ) is a global weak solution 
of (12. 8h in if and only if is in W^^°(Q) and satisfies equation (12. 8p and the 
Rankine-Hugoniot condition on 5": 

[p(|Dvf ,^LV-i% = 0, (2.15) 

where [F]s is defined by 

[F(t,ri)} s := F&rj)]^- F&rj)^ for (^)e5. 
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Note that the condition <p G W l( £°(£l) requires 

[tp] s = 0. (2.16) 

In Figure [TJ by (|2.12p and (|2.16p . the pseudo-potentials ipo and <pj\f below So and 
Sj\f are respectively in the form of 

fo = ~7;(f + V 2 ) + u o£ + voV + k o, 

\ (2-17) 

<PAf = ~2^ 2 + + u ^ + V -WV + k M 

for constants uo,vo,uj\fi W> ko, and kj^. Then it follows from ()2.9[) and (|2.17p that 
the corresponding densities po and pj\f below So and Stf are constants, respectively. 

Given > 1, we obtain (uo,vq) and by using the shock polar curve for 
steady potential flow in Figure [2 In Figure [21 let 6 son i c be the wedge angle such 
that the line ^ = tan 6 SO nic intersects with the shock polar curve at a point on the 
circle of radius 1. For a wedge angle 9 W 6 (0,9 son i C ), the line v = ut&n6 w and the 

shock polar curve intersect at a point (uo,vo) with \J u + v a > 1 and uo < Uoo. 
The intersection (uo,vo) indicates the velocity for steady potential flow behind an 
oblique shock So attached to the tip of the wedge angle 9 W . Since the strength 
of the shock So is relatively weak compared to the other shock given by the other 
intersection point on the shock polar curve, So is called a weak shock. We also 
note that the states on both sides of So are supersonic, and such states (uo,vq) 
smoothly depend on Uqo and 9 W . Once (uo,vo) is determined, then, by (|2,13p . the 



V 

t 




Figure 2. Shock polar 



corresponding density po is given by 

Po'' = 1 + " «o - vl). (2.18) 

To find a curved shock S connecting the oblique shock So with the normal shock 
Sj\f in Figure [lj it is more convenient to use a coordinate system in which the 
boundary of wedge V we d ge is the same for all 9 W 6 (0,9 son i c ). For that reason, we 
introduce a new coordinate system (£',?/) defined by 




(2.19) 
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In other words, (£', r/) is obtained by rotating (£, 77) and translating the self-similar 
plane. Denoting as (£,77) for v?oo and (po defined by (|2.12p and (|2.17p are 

expressed as 



foo = ~\^ 2 + V 2 ) - Woo sin6 w + k c 
fo = ~t:(£ 2 + V 2 ) + ubi + k o 



(2.20) 



with 



UO = \J Uq + Vq — Mqo cos 9 W (2-21) 

for constants and ko different from (|2.12[) and (|2.17|) . For simplicity, we will 
write uq as uq hereafter. Without loss of generality, we choose koo = 0. Set 



^00 := sin I 



(2.22) 



and let /3 be the angle between the oblique shock So and wedge T we d ge - We use the 
parameters (i>oo>/3) instead of (uoo,9 w ) to compute the normal shock SV and the 
state behind it. 




Figure 3. Weak shock solutions on the self-similar plane 



By (j2.16j) . the definition of (poo given by (12.2QH with koo = 0, and by the slip 
boundary condition (|2.14j) . if we let Stf = {rj = r/jv"} for a constant rjj\f, then we 
have 

-l(f + V 2 )-VooW- (2.23) 



By (|2.15p and (j2. 13j) . the density pj\f and constant rjjy satisfy 

Voo 



W = 

7-1 



PN - 1 
- 1 1 



7-1 



^"00 + W^oo- 



(2.24) 
(2.25) 



We note that < r\^[ < cj^f for cj\f = pfy . 

Fix j3 E (0, f ), consider an oblique shock So of the angle j3 from the £— axis, 
and let tpo be the corresponding pseudo-potential below So- Let (£o,Vo) De the 
intersection Pi of So and the sonic circle B co (uo,0) for the state in OP (see Figure 
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1). There are two intersection points, and we denote Pi the point with the smaller 
value of the £— coordinate (see Figure [3]). Given Vqo > and 7 > 1, rjo satisfies 

^-<0 for all p G (0, |). (2.26) 

Therefore, the set I Voo := {(3 G (0, ^) : < rjo(P) < 00} is connected and, further- 
more, there exists (3 son i c G (0, ^] satisfying 

^«oo — (0 1 ft sonic) j 

and such j3 SO nic depends smoothly on and 7. 

Given 7 > 1, we note that sonic on the shock polar curve in Figure [2] depends 
on Uoo, while f3 son ic depends on v^, so that we write as 6 , sonic (u 00 ) and f3 S onic(voc), 
respectively. Define two parameter sets *Pi and ^2 by 

?Pl := (J {^00} x (0, 

"sonic (Uoo)), 

Uoo>1 (2.27) 

:= |J {«oo} X (O^gdic^oo)). 

Voa>0 

Then we have the following lemma. 

Lemma 2.1. For any 7 > 1, there exists a one-to-one and onto correspondence 
between *Pi and ^2 ■ 

Remark 2.2 (Normal shock: /3 = 0). For fixed 7 > 1 and > 0, the state of 
(3 = is a normal shock, and this corresponds to the state of (itoo,0u>) = (00, 0). 
Even though {3 = is nonphysical, we will put the case of ft = in our consideration 
as it is useful in applying the Leray-Schauder degree theorem to establish Theorem 
[H stated in Section [3l 

For a fixed (voo,(3) G if we can prove the existence of a curved shock r^o,^ 
to connect So with S^ and the existence of corresponding pseudo-potential below 
F shock, then, by Lemma [2~T| the existence of a curved shock and the pseudo-potential 
below the shock for (uoo,9 w ) G *Pi corresponding to the fixed (voo,/3) automatically 
follows. By taking (voo,(3) as a parameter instead of («oo >#«))> it is advantageous 
that the boundary of a wedge is always on the £-axis. As (voo,f3) is the parameter 
of our problem, we write ipoo in (|2,20p as 

Voo = -\(? +r?)-v 00 r ] . (2.28) 

Since the states below So and S>/ are given by ipo and ipjij, in order to solve 
Problem O it suffices to solve the following problem. 

Problem 3 (Free Boundary Problem). Find a curved shock T s h oc k & n d a function 
(p defined in the region f2, enclosed by T shock, ^fonic^sonio an< ^ iv = 0}; su °h that 
ip satisfies 

(i) Equation (|2.8p in Q; 

(ii) (p = (foo, pDip ■ u s = Dtfoo ■ v s on T shock ; 

(iii) if = ipp, Dip = Dipp on Tg nic U T{f onic for <p p := max(ipo, <p/f); 

(iv) d„ip = on T wedge , 



8 



MYOUNGJEAN BAE, GUI-QIANG CHEN, AND MIKHAIL FELDMAN 



where v s is the interior unit normal on T s h oc k- 

Let ip be a solution of Problem [3] with a shock T shock- Moreover, assume that 
if G C 1 (Q), and T s h oc k is a C 1 -curve up to its endpoints. To obtain a solution of 
Problem [2] from ip, we divide the half-plane {n > 0} into four separate regions. Let 
Q.E be the unbounded domain below the curve So U T s h oc k U <SV in {?7 > 0} (see 
Figure [3]). In fig, let f2e> be the bounded open domain enclosed by So,^® on i c , and 
{77 = 0}. We set Clj\f := Qe\ (^e> U Define a function 92* in {77 > 0} by 



Poo 


in {7/ >o}\n E , 




in 




inr ■ lj mi ■ 

m 1 some u J L u 1 sonic' 




in Q_\f. 



(2.29) 



By (12.160 and (iii) of Problem [3j 99* is continuous in {77 > 0} and is C 1 in Qe- In- 
particular, 99* is C 1 across r^UT^. Moreover, using (i)-(iii) of Problem [3l we 
obtain that p* is a weak solution of equation (|2.8j) in {7/ > 0}. Applying the inverse 
coordinate transformation of (|2.19j) to p*, we obtain a solution of Problem [2 



3. Results 

In this section we first present the known result, Elling-Liu's theorem, in [13J, 
and then we describe our new results for Problem 3, the Prandtl-Meyer problem. 

3.1. Known Result. Elling and Liu in [13] have established the following theorem. 

Elling-Liu's Theorem (Theorem 1 in [13]). Given (voo,/3) £ ^2> let L be the line 
segment connecting Pi with P2 in FigureO If L does not intersect with Pi(0, —v^), 
then there exists a weak solution to Problem [2] with structure (12.290 . 

If Voo > 1, since Pi(0, — Voo) C {n < 0}, the assumption in the theorem above 
holds true. On the other hand, for < 1, there exists a set of parameters (voo, (3) S 
^2 for which the assumption no long holds. This can be shown as follows: 

Fix 7 > 1 and Voo > 0. Let be the line with positive slope such that 
is tangent to the sonic circle Pi(0, — t>oo) of the incoming state and that it passes 
through the point P2 = rjj^f) = T^ onic n Sj\f. For some Voo > 0, there may be 
two such tangent lines. In that case, we choose the one with the smaller slope. We 
also fix /3 € (0, (3 S onic(voo))- Let Lq be the line segment connecting Pi := (£e>, r/o) = 
fg~D So with P 2 (see Figure H}. 

Let tan^o and tan be the slopes of the lines Lo and Loo, respectively. If 
®o > $00 > then the line Lo intersects with the sonic circle Pi(0, — 7>oo)- On the other 
hand, if 80 < 0oo, then Lo has no intersection with Pi(0, — Woo)- Define 

P(^) := tan6»o -tan6»oo. (3.1) 

Then, we have the following proposition. 
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Figure 4. Top: o < 9oo] Bottom: O > #00 

Proposition 3.1. For any given 7 > 1, there exists v * > depending on 7 so that, 
wherever Voo E (0, i>*) ; there exists f3 = /3(t>oo) G (0, f3 S onic(voo)) satisfying 

F((3)<0 forpe[0j], 

F((3)>0 for f3e 0,/3 sonic ). ( °' _) 

According to Proposition 13.14 there exists a subset of ^2 for which (131 Theorem 
1] does not apply. 

3.2. New Results. Our new results provide not only the existence of a global 
solution to Problem [3] for all (voo,/3) in but also the higher regularity of the 
solution in pseudo-subsonic regions. Therefore, we achieve the existence of a self- 
similar weak solution with higher regularity to Problem [1] for all (uoo,6 w ) in 

For any fixed (voo,f3) E ^2, the oblique shock So of the slope tan/3, the normal 
shock SV) and the corresponding pseudo-potentials ipo and (p/s below So and Sj\f 
are uniquely computed. Let (£(/?), 0) be the ^-intercept of So- Then we have 

So C {(£, f o (0) ■ fo(0 = tan/?(£ - G M}. 

Denoting Pi = (£e>,r?e>) and P2 = (£aT> w) m Figure O we have 

So = {(£, r?) : £(/?) < 6 < Co, r/ = /©(£)}, SV = {(£, W) : ? > Cat}- 

Theorem 1. For any (/wen 7 > 1 and (woo>/3) G *$2; i/iere exists a global weak 
solution ip of Problem [3] satisfying the following properties: 

(i) There exists a shock curve T s h oc k with endpoints P\ = (£0,770) an d P2 = 
w) such that 

- ^ shock satisfies T shock C (R% \ B^O, -Uqo)) n {£e> < £ < 6\f}; 

- rs^ocfc is C 3 m ife relative interior: For any P E ^ s hock\{Pi, P2}, there 
exist r > and f G C 3 (1R) suc/i i/iai 

r s/locfe n B r (p) = {(£,//) : t? = /(£)} n b p (j>) ; 
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- The curve Sq U ^ shock U Sj*/ is , including at the points P\ and Pi', 

- ^shock^onic^Znic and T wedge := {(£,0) : u Q - c a < £ < c^} do not 
have common points except at their end points. Thus, T s h ck U ^onic U 
^sonic U ^ wedge is a closed curve without self-intersections. Denote by Q. 
the open and bounded domain enclosed by this closed curve. 

(ii) <p satisfies p G C 3 (H \ (Tf onic U I^J) D C\Tl). 

(iii) Equation QZJft is strictly elliptic in fi \ (r^ nic U rf onic ). 

(iv) max(yo,^) < ^ < v?^ m f2. 

(v) <9_|(</?oo -<p) < and <9e So (^oo - y?) < in /or 

If 99 is a weak solution satisfying properties (i)-(v) of Theorem [TJ then </? and 
^ shock satisfy additional regularity properties. 

Theorem 2. Given 7 > 1 and (voo,f3) E ^?2; tet <p be a weak solution of Problem Q] 
satisfying properties (i) -(v) of Theorem Q} Then the following properties hold: 

(i) The curve So U r s /j OC /% U SV zs C 2,a for any a G [0, i), including at the 
points Pi and Pi. Moreover, Y "shock is C°° in its relative interior. 

(ii) <p g C°°(n \ (r? onic u r^J) n C 1 - 1 ^). 

(iii) For a constant a > and a set P given by 

v = {(€>V) '■ m ^(Po(^,v),PAf(^v)) < cpoo(£,v),V > °}> 
de/me 

P^ = Pn{(£,n) :dist((C,r?),r^ nic ) < a} n B CM (0,0), 
V° =Vn{{t;,r } ) :dist((e,r?),rf onic ) <ff}nB C£) (« Ol 0) 

/or CjV = PyJ and cq = Pq ■ Then, for any a G (0, 1) and any given 
(Co,??o) G (r^ nic U T^ nic ) \ {Pi,P 2 }, there exists K < 00 depending only on 
7,foo, e 0) a, [Mle*.i(nn(z>fux>^))> anrf d = dist ((£o, Vo), T sho ck) so that there 
holds 

^k,a,nnB d/2 (So,vo)n(v? Q/2 uv° /2 ) ^ K ' ( 3 ' 3 ) 

(iv) For any (£„, %) G (lf onic U T^ nic ) \ {Pi, ft}, 

lim (A-rV 9 - -Drr max^jv", 990)) = — — -, (3.4) 

(«,j))-»(eo^o) 7 + 1 

ui/iere r = y/£ 2 +n 2 near r^ nic , and r = y^(| - uq) 2 + n 2 near T^ c . 

(v) The limits lim L> 2 </? and lim D 2 gj do not exist. 

Remark 3.2. Assertion (iii) 0/ Theorem [2] says tfoat ip in Q is C 2,a up to the sonic 
arcs T® onic and Y{f onic away from the points Pi and P2 . Assertion (ii) of Theorem 
[2] combined with (iii) 0/ Problem [3] imply that the function p* in (|2.29p is C ' 
across the sonic arcs T^ nic andVf onic . Assertion (iv) 0/ Theorem [2] implies that the 
function is not C 2 across the sonic arcs r^ nic and T® onic since there is a jump 
of second derivative of <p, which implies that the C 1 ' -regularity is optimal. 
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Remark 3.3 (Weak solutions of Problem [3] for /3 = 0). For (3 = 0, ip = tp^f is the 
unique weak solution of Problem [3] satisfying properties (i)-(v) o/TheoremQ], where 
cpj\f is defined by (|2.23p . 

The following theorem easily follows from Theorem [IJ Lemma |2.1| and the argu- 
ment after (pl29j) . 

Theorem 3. For any given 7 > 1 and any (uoo,9) G ^Pi, Problem [2] has a global 
weak solution p* of the structure as in Figure Q] with ip* being continuous in A and 
C 1 across T° onic U T^ onic . 



4. Overview of the proof of Theorem [T] 

The key to establish Theorem [T] is in the monotonicity properties of poo — ip for a 
weak solution ip of Problem [3l It implies that the shock is Lipschitz graph in a cone 
of directions and thus fixes geometry of the problem, among other consequences. 
Another key property is following. Fix 7 > 1, t>oo > 0, and /3* G (0, (3 son ic{v 00)) ■ 
Then there exists a constant C > depending only on 7,1*00, an d /3* such that, for 
any f3 G (0, /?*], a corresponding weak solution 92 with properties (i)-(v) in Theorem 
CD satisfies 

dist(T shock , BxiO, -Voo)) > — > 0. 

This inequality plays an essential role to achieve the ellipticity of equation (|2.8|) in 
f2. Once the ellipticity is achieved, then we obtain various apriori estimates of <p so 
that we can employ a degree theory to establish the existence of a weak solution for 
all (3 G (0, /3*]. Since the choice of /3* is arbitrary in (0, (3 SO nic(v 00)), the existence of 
a weak solution for any (3 G (0, (3 son i C (voo)) can be established. 

For the rest of this paper, we outline the proof of Theorem [TJ 

4.1. Monotonicity Property and Its Consequences. Given Voo > and (3 G 

(0,/?*], set e~s := (cos /3, sin/3). Then e~s ls a un it tangent to So and satisfies 
e So ■ i > 0. 

Lemma 4.1. For a fixed (voo>/3) G if p is a weak solution with properties 
(i)-(v) in Theorem [T], then it satisfies 

d eSo (p oo -p)<0 inU\f°~, d^(p oo -p)>0 ron\r£~, 
£^(9? - <px) < 0, £^(9? - px) < in J7, 
9e So (p - P>o) > 0, d v (p - po) < in fl 

For /3 G (0, /?*], and two unit vectors e»s and — £ in M 2 , define 

Cone(<55 ,-£) := {aie 5o +«2(-|) : ai,a 2 > 0}. 

Note that the vectors es a and — £ are not parallel if (3 G (0, /?*], thus, Cone(e5 , — £) 
has non-empty interior. Let Cone (es o , —£) be the interior of Cone(es , — £). 
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Remark 4.2. For (3 = 0, we have es = £. Then, for (3 = 0, we define Cone(is , — £) 
by the upper half-plane, that is, 

Cone(e So ,-£) := {u G M 2 : u ■ f] > 0}, 

and let Cone (es o , — £) 6e i/ie interior of the upper-half plane {it G M 2 : u ■ fj > 
0}. Note that this is consistent with the definition for j3 > in the sense that 
(Cone (e,s o , — £))|/3 — > (Cone°(es i c , (/3 = 0),— £))\p =0 as ft — > 0+, where the conver- 
gence is locally in the Hausdorff metric. 

Hereafter, we assume that <p is a weak solution satisfying properties (i)-(v) in 
Theorem [H unless otherwise specified. By Lemma 14.11 we have 

d e (<Poo — (f) < in SI for all e G Cone (es o , — £), 

<9 e (</? — CjCjv) > in $7 for all e G Cone(— £, —77), (4.1) 

d e (ip — tpo) > in O for all e G Cone(es , —77). 

By (|2.16p and the first inequality in (|4.ip . if ip is a weak solution with properties 
(i)-(v) in Theorem[T]for (voo, /3) G there exists a function 77 = f s hock(€) satisfying 

(i) r s/l0cft = {((,!]) : £ G (tp^tpj,?? = /s/iocjfc(0}> where ^ is the ^-coordinate 
of the point -Pj for j = 1, 2; 

(ii) fshock(0 is strictly monotone in £ G (£o,£/v): 

£/kh*(0>0 for ail £€(&>, (4-2) 

(iii) there exists a constant C\ depending only on 7, v^, and /3* such that 

sup |/^ ocfc (OI<Ci (4-3) 

From (|4.2p . it easily follows that 

inf dist(T shock ,T wedge ) > inf rj Pl > 0. 

/3G[0,/3*] [OA] 

Moreover, f|4.2j) implies that the region Q enclosed by r sftocfc , r£ nic , T^ nic , and 
r«ie<ige is uniformly bounded for all (3 G (0, /?*]. In fact, we have 

O C {(£, 77) : ue, - co < £ < c/v, < 77 < t/aa}. 

From this, we obtain a constant C > depending only on 7, foo an d /3* such that, 
if </? is a weak solution with properties (i)-(v) in Theorem Q] for /3 G (0, /?*], then 

C B o (0), 

supM < C, |M| c o,i(n) < c » 
n y ' (4.4) 

2 _l_ 

( — — r) 7 " 1 </5<CinS7, 1 < p < C on Y shock . 
7 + 1 

The importance of the last two inequalities in (|4.ip will be mentioned in Section 
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4.2. Uniform lower bound for the distance between the shock and the 
sonic circle of the upstream state. We can establish the following crucial propo- 
sition. 

Proposition 4.3. Given 7 > 1 and Voo > 0, if tp is a weak solution for j3 G (0, /3*] 
with properties (i)-(v) in Theorem [TJ then 

dist(r shocfe , B^O, -««,)) > ^ > 0. (4.5) 

Proposition 14.31 is essential to obtain uniform ellipticity constants of equation 
(|2.8p . More precisely, we first write equation (I2.8P with p given by (|2.13p and 

Poo = 1 as 

divA(D<p, <p, e, v) + V, £,V) = ( 4 -6) 

for p = (pi,p2) G M 2 , z G M, and (£,77) G M 2 , where 
A(p, z, £, V ) = A(p, z) := p(\p\ 2 , z)p, B(p, z, £, 7?) = B(p, z) := 2 /9 (|p| 2 , *). (4.7) 

Then, using (|4.5p . we can find a constant C such that, if ip is a weak solution for 
/? G (0, /?*] and a set J7 C f2 satisfies dist(?7, F® onic U r^ nic ) > d > 0, we have 

d 2 

^\k\ 2 < A;.(D^,r,),^,r,)) KiKj < C\k\ 2 (4.8) 

for (£,77) G 17 and any k = (^1,^2) G M. 2 . The important part is that such a 
constant C can be chosen depending only on 7,^00, and /?*, but independent of 99 
and therefore we obtain the uniform estimates of weak solutions of Problem [3] for 
/3 €((),&]. 

4.3. Uniform Estimates of Global Weak Solutions. We fix v^o > 0, 7 > 1, 

and /3* G (0, /3 S onic(voo)) ■ Thanks to (I4.8p . we can achieve uniform estimates of 
weak solutions ip in Q for /3 G (0, /3*]. Because of the degeneracy of the ellipticity of 
equation (|2.8p on r^ nic U r^ nic , we need to consider two cases: 

Case 1: Away from Tf onic U r^ nic , where equation (]2.8p is uniformly elliptic; 
Case 2: Near r^ nic U T^f onic , where the ellipticity degenerates. 

4.3.1. C 2,a -estimates away from r^ nic ur^ nic . By (|4.8|) . equation (|2.8|) is uniformly 
elliptic away from the sonic arcs Tf onic U r^ nic . Then we can show that, for any 
a G (0, 1) and r > 0, there exists a constant C > depending only on 7, Voo, a, 
and r such that, if is a weak solution with properties (i)-(v) in Theorem [T] for 
P G (0, /?*], then we have 

(i) for any B 2r {P) C O, 

H^llc 2 ' Q (B^P)) - C '' 

(ii) for P G r w , e and (B 2r (P) nfi)n (T^ U I^J = 0, 
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For uniform estimates of shock curves T^odfe and weak solutions <p on r^^, we use 
a partial hodograph transform. For any f3 G (0, /3 son i c ), we have 

fj G Cone (e5 o ,-|). 

Then, by the first inequality in (|4.ip . we can find constants 5 and a > depending 
only on 7,Uoo ; and /3* such that, if ip is a weak solution for /3 G (0,/3*], then 

^(<^oo - v>) < in -^(r^ocfc) n n (4.9) 

for N a (Pshock) = {(£,V) G ^ 2 : dist((£,r]),T shock ) < a}. We introduce a coordinate 
system (£',?/) with £' = £ and a function v(£',r]') such that u satisfies 

u(£,(poo-p)(M) =r?. (4.10) 

By (|4.9|) . such a coordinate system is well defined. In the new coordinates, 

we use equation f|2.8j) in the set J\f a (F shock) and the Rankine-Hugoniot condition: 

pDip ■ v s = Dip^ ■ v s on T shock (4.11) 

to obtain an equation and a boundary condition for v. Then we obtain an ellip- 
tic equation and an oblique boundary condition for v in the (£', r/)-coordinates. 
By (|2.16p . Ts/jocfc becomes a fixed flat boundary; then we obtain the uniform esti- 
mates of v for j3 G [0, /?*]. By f)4. lOj) . we obtain T s h oc k = {r] = i>(£,0)}; thus the 
uniform estimates of v imply that, for any k G N and d > 0, there are constants 
s,Ck(d), and Ck{d) depending on 7, t>oo>/3*, and d such that, if P(Cp,i]p) £ r s h oc fc 
and dist(P, Tf omc U T^ mc ) > d, then 

|£ fc /.h«*(&>)l < Cfe(d), I^LH ^ ^) 011 ^( p ) n °> ( 4 - 12 ) 

where Ck(d) and Ck(d) depend additionally on 

4.3.2. Weighted C 2 ' a -estimates near T® onic U r^ nic . Near the sonic arcs r^ nic U 
r^Lv c , it is convenient to use the coordinates in which the sonic arcs are flattened. 
For that reason, we introduce the (x, y)-coordinates as follows: Choose e > small. 

• For (£,77) G A4(r^ mc ) := {(£,77) : dist((£, 7?), T* mc ) < e}, we set 

x := cj\f — r, y := 9 

for the polar coordinates (r, 9) centered at (0,0); 

• For (£,7?) G A4(r? omc ) := {(£,7?) : dist((£, 77), T^J < e}, we set 

x := Co — r, y := tt — 9 

for the polar coordinates (r, 9) centered at (uo,0). 
From the definition of (x, y) above, it is easy to see that 

(finA4(r)) c {x > o,y > 0}, r = Hn{x = o}, 
({ ?? = o}nA4(r)) = nn{ y = o} 

for T = Tf onic or T^ onic . For and tpj\f defined by (j2.20R and (I2.23|) . respectively, 
set 

i/):=<p — max(^o, ip^). 
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In fact, choosing e sufficiently small, we have 

jipo m A4(rf omc )nft, 

t^W m A/ £ (F s v onic ) n 0. 

This implies that 

rt p = ip - (po i n A4(rf onic ) no, ij) = ip — <ptf in ■Af e (r££ nic ) n a 

Since we seek a weak solution Lp of Problem [3] to be C 1 across r^ nic U Tf onic , ip 
satisfies 

^(0,y) = 0. (4.13) 

Using the definition of ip, we can rewrite equation (|2.8p as an equation for ip in the 
(x, y) -coordinates near Tf onic or r^ nic as follows: 

An(Dip, ip, x)ip xx + 2A 12 (Dip, ip, x)ip xy + A 2 2(Dip, ip, x)ip yy + A(Dip, ip, x) = 0, 

where (Aij, A)(p, z, x) are smooth with respect to (p, z, x). Then (|4.8p implies 

x I ,2 ^ A xl (Dip,ip,x) 2 A 12 (Dip,ijj,x) 2 1 2 

A ft < — KiK 2 + ^22 (-DV'i V,X)Kn < T\ K \ ( 4 -14) 

x Jx A 



for a constant A > 0, where A depends only on Vqo, 7, and /3*. Using the expressions 
of the coefficients Aij(p,z,x) and the fact that an ellipticity constant in (|4.8j) is 
proportional to the distance to the sonic arcs T® onic U r^ nic , we obtain a constant 
5 > such that, for any /3 G (0,/3*], a corresponding ?/> satisfies 

^ < 7X^^ in Wfonic U T^ mc ) n 0, (4.15) 
1 + 7 

and such <5 > depends only on 7,^00, and /?*. Also, by the last two inequalities in 
()4.ip and the definition of the (x, y)-coordinates, one can easily check that ip x > 
in A4(r^ nic U T^ nic ) n fl Combining this with (|4.15p . we obtain 

|^| < Cx in o n A4(r^ nic u rf onic ). (4.16) 

Combining (|4.16p with (|4.13p . we obtain 

Mx,y)\ <Cx 2 in ^(iifcUlljna (4.17) 

It is important to note that the constant C in (|4.16|) and ()4.17|) depend only on 
7, t>oo, and /3*. Then the uniform weighted C 2 ' Q -estimates of ip for j3 G (0,/?*] are 
obtained by the local scaling of ip and the covering argument. 

For P E A4(r? onic U T^ mc ) n (fi \ (rf onic U r^ nic )), set P := (*o,J/b) in the 
(x, y)-coordinates and define 

^o,yo) (5 T) = J_^ (xo + djgr yo + with d = ^0 

er 1U 

in Qj*' '^ := {(S,T) G (-1, l) a : (x + dS,y + v^T) G fi}. For the estimates of 

H^O^o.ifo) || — — ^ ; we need to consider the case of Po E and the case of 

C 2 ' Q (Q 1(/2 ) 

Po ^ Ts^ocfc, separately. For Po G r^cfe, we use the first inequality in (14. ip to show 
that 

d y (tpoc -<p)>5 near rf omc or T^ onic for some <J > 0. 
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Then, as in (|4.10p . we can introduce a coordinate system (x',y') and a function 
w(x',y') satisfying 

w{x, - tp)(x,y)) = y, 

and we write an equation and a boundary condition for w by using equation (|2.8|) 
and the boundary condition (|4.1ip . And, we take the same procedure to the function 
ip oo — max(ipo, (pj\f) to introduce a function wq and its equation. Then we take the 
difference of the equations for w and u>o to obtain an equation for w — Wq and use 
it to estimate w — Wq. This procedure will provide the weighted C 2 ' a -estimates of 
Y shock near P and ||^o,,o) ||^ . 

For Pq G Q U T we d ge , it is easier to obtain the estimates of \\ifj( xo >y°) || (x , vo ) ■ 

c 2 ' a (Qi /2 ' ) 

Then, the supremum of \\^.vo)\\^ a ^—^ ) over (x , y ) G M^^T^) n 
is uniformly bounded for any (3 G (0, /?*]. We define 

ll^llSicDnn == f -p o U {x °' y \ 2 , a{ ^ ) 

for r = Tf onic U T^ nic . We note that the uniform estimate of IMIj££v e(r ) nn auto " 
matically provides the uniform C 1,1 -estimates of ip in Af E (T) n for /3 G (0, /3*]. 

4.4. Outline of the proof of Theorem [lj, To establish Theorem [TJ that is, the 
existence of a weak solution of Problem [3] for all (vqo, 0) G ^2, we apply an iteration 
procedure. For that purpose, we define an iteration set K, and an iteration map T 
on /C so that a fixed point of T yields a solution of Problem [3l Then, we apply the 
Leray-Schauder degree theorem to show the existence of fixed points. 

Given 7 > 1 and (voo,P) G let <p be a solution of Problem [3] with a shock 
r shock satisfying (i)-(v) of Theorem [TJ By (ii) of Problem El T s h oc k is determined 
by the solution if. Moreover, by the monotonicity property (|4.ip . T s h oc k is a graph 
in ^-direction, so we can introduce an invertible mapping & which maps onto 
the rectangular domain Q lter ~ (—1,1) x (0,1). Then the function U = pS" 1 
is well defined on Q lter regardless of fi. In other words, any solution of Problem [3] 
satisfying (i)-(v) of Theorem [TJ can be re-defined as a function on Q lter . This allows 
us to perform iteration in a set of functions defined in Q lter to prove the existence 
of a solution of Problem [3l 

We define the iteration set /C as a product of a parameter set and a set of functions 
defined in Q lter . Given 7 > 1 and Voo > 0, we fix /3* G (0, (3 SO nic(voo)) and set 

K := U moM {f3} x /C(/3), 

where /C(/3) is a set of functions defined in Q lter . We define /C(/3) such that, if 
u G JC(/3), there exist a corresponding pseudo-subsonic region fiW an d a function 
defined in fiW, which satisfies several properties including the monotonicity 
properties (|4.ip . so that it can be a candidate of solution of Problem [3] satisfying 
(i)-(v) of Theorem [TJ 

Once the iteration set /C is defined, the iteration map J- is defined as follows: 
Given ((3,u) G /C, define the corresponding pseudo-subsonic domain Q( u \ enclosed 
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b y ^onic^Znia and ^ wedge, and an approximate shock T^ ock , and set up a bound- 
ary value problem for an elliptic equation whose ellipticity degenerates on T{f onic U 
^?onic- Let (p be the solution of the boundary value problem in Q( u > . The iteration 
set /C will be defined so that such dp exists. Then we can define a function u in Q lter 
from dp. Then the iteration map T is defined by dF(/3, u) = u. The boundary value 
problem for dp is set up so that, if u = u, then dp satisfies equation (|2.8j) . with an 
ellipticity cutoff in a small neighborhood of F® onic U T^ nic , the Rankine-Hugoniot 
conditions (|2. 15[) — f)2. 16[) on the shock r s ^ oc fc, and the boundary conditions stated in 
(iii)-(iv) of Problem[3l More specifically, note that two conditions in (iii) of Problem 
[3] are specified on the sonic arcs. Since the sonic arcs are fixed boundaries, this looks 
like an overdeterminacy. Thus, for the iteration problem, we use only the Dirichlet 
condition 

dp = ipp = max( Wl tpjf) on rf onic U T^ mc , 

and then we prove that the condition Dip = D<pp on Tf onic U T^f onic also holds. In 
this proof, we use the elliptic degeneracy of the equation in fiW near r2, nic ur^ nic . 
Technically, this follows from the estimates of ip = <p — pp in norms (I4.18P near 
Y° ■ U ■ . 

sonic sonic 

If (/?,«*) G K. satisfies F{f3,u*) = u*, then we call u* a fixed point. In order 
for a fixed point to provide a solution of Problem [3j we need to show that p>* 
satisfies equation (|2.8p without the ellipticity cutoff. Moreover, since we intend to 
apply the Leray-Schauder degree theorem to establish Theorem [TJ we need to prove 
the following properties: 

(i) If u* G /C(/3) is a fixed point of T for some j3 G [0, /3*], then the corresponding 

satisfies equation (]2.8h without the ellipticity cutoff; 

(ii) For any j3 G [0, the map T(f3, •) is compact, and T is continuous; 

(iii) The iteration set fC is open; 

(iv) For any /3 G [0,/3*], there is no fixed point of T on the boundary of the 
iteration set IC(f3). 

In proving all the properties above, the apriori estimates in Sections 4.1 and 4.2 
play an essential role. Once (i)-(iv) are verified, then deg(J r (/3, •) — Id,IC(f3),0) is a 
constant for all (3 G [0, 

To complete the proof of Theorem[[J we show that deg(J r (0, -) — Id, /C(0), 0) ^ in 
two steps. First, we prove that ^(0, •) has the unique fixed point u ( normal ) i n /C(0). In 
fact, we have u ( norma ^ = 0. Then we can conclude that deg(J c '(0, -)-Id, K(0), 0) ^ 
by showing that D u J r (0,u ( - norma ^) — I has the trivial kernel. 

As mentioned earlier, once we establish Theorem [H that is, the existence of global 
weak solutions of Problem El we conclude the existence of a global self-similar weak 
shock solution for any admissible wedge angle 9 W by the one-to-one correspondence 
between the parameter sets and 

For the detailed proof, see Bae-Chen-Feldman [2]. 
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